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The earlier paper of this same title demonstrated a collapse of relaxation data of fragile supercooled 
glass forming liquids [J. Phys. Chem. B 113, 5563-5567 (2009)]. For temperature T below that 
of the onset to supercooled behavior, T , the logarithm of structural relaxation time, logr, is given 
by the parabolic form log(r/r ) = J 2 (l/T— 1/T D ) 2 , where J and r Q are temperature independent. 
This paper presents further applications of this formula. In particular, it is shown that the effects 
of attractive forces in numerical simulation of glass forming liquids can be logically organized in 
terms of J and T . Further, analysis of experimental data for several systems suggests that J and 
T are material properties. In contrast, values of similar parameters for other fitting formulas are 
shown to depend not only upon the material but also upon the range of data used in fitting these 
formulas. Expressions demonstrated to fail in this way include the Vogel-Fulcher-Tammann formula, 
a double-exponential formula, and a fractional exponential formula. 



I. INTRODUCTION 

Transport properties of supercooled glass forming liq- 
uids are strong functions of temperature, increasing by 
typically ten or more orders of magnitude while abso- 
lute temperature is reduced by only ten or twenty per- 
cent [TJ |2] • Most often this growth with respect to low- 
ering temperature is super- Arrhenius and materials that 
behave in this way are called fragile. Seemingly varied 
behaviors of fragile glass formers are catalogued in terms 
of a property known as fragility [3]. This property mea- 
sures the relative rate of change of a transport property 
as temperature is reduced near the glass transition tem- 
perature T g - the point at which the material falls out of 
equilibrium. Fragility varies significantly from one mate- 
rial to another [4]. Yet, despite this variability, we have 
found that all such transport data can be collapsed to a 
non-singular function of temperature T jS] , 

log (r/r ) = J 2 (l-^) , T<T , (1) 

where r stands for the transport property, such as struc- 
tural relaxation time or viscosity (which may also be de- 
noted as 77), and r D refers to that same property at the 
onset temperature, T . 

The onset temperature marks the crossover from 
normal liquid behavior to supercooled liquid behavior. 
Above that temperature, transport is unremarkable, in- 
deed nearly temperature independent [SHI]. It is only 
below that temperature where most equilibrium liquid 
transport properties depend sensitively upon tempera- 
ture. It is generally accepted that this crossover is due to 
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the appearance of local rigidity or dynamical constraints, 
and this feature makes molecular motion rare and acti- 
vated |H [10]. It is similarly accepted that correlations be- 
tween motions in different regions of space cause activa- 
tion barriers to grow with growing length scale, and this 
growth produces super-Arrhenius behavior for T < T Q 
(for example, Refs. 2, 9 12 ). But the mechanism and 
nature of this growth remain unresolved. 

To reach a resolution, it would be helpful to know if 
there is a precise and general form of temperature depen- 
dence for supercooled liquid transport. From this per- 
spective, it seems significant that throughout the range 
T g < T < T , where transport properties vary by many 
orders of magnitude, all equilibrium transport of frag- 
ile glass formers is described by Eq. [I] This formula 
is an exact result for a class of dynamical models [T3T - 
I17j . Can other formulas based on other models prove 
equally satisfactory? We address this question in this 
paper by comparing Eq. [TJ with a few other proposed re- 
lationships. The alternatives are the commonly employed 
Vogel-Fulcher-Tammann (VFT) form [IS], the recently 
advocated double exponential form |19j , and a variant of 
Eq. [TJwith a fractional exponent less than 2 [20]. Each 
of these, like Eq. [TJ have three parameters. (We exclude 
formulas with four [21] . or more, parameters from our 
consideration.) We show that the parameters of Eq. [TJ 
are material properties in the sense that their values for 
a given system are reasonably invariant to the number 
and range of experimental data points. The parameters 
associated with the other considered forms do not have 
this quality. 

The fact that data can be collapsed with Eq. [TJ using 
parameters that are stable with respect to the range of 
data implies that this relationship can be used to pre- 
dict results of measurements not yet performed. It also 
provides a means to sensibly organize data. Xu, Liu and 
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TABLE I: Comparison between VFT and parabolic fits for supercooled liquid OTP 





VFT" 


Parabolic 11 


Fitted temperature range 


l°g^vFT/ Poise 


A/K- 1 


TVft/K 


Error 6 


J/K- 1 


T /K 


log r\ j Poise 


Error 6 


239 K - 350 K 


-8.98 


1.03 • 10 3 


191 


0.20 


2.79 • 10 3 


352 


-1.88 


0.064 


267 K - 350 K 


-4.30 


254 


238 


51 


2.82 • 10 3 


351 


-1.86 


0.087 


239 K - 265 K 


-8.38 


999 


191 


0.29 


2.82 • 10 3 


349 


-1.60 


0.13 



a listed fitting parameters are determined by minimizing the mean square deviation between the fitting function and experiment for the 
specified range of data. 

6 Root-mean-square deviation between fitted function and experiment over entire set of supercooled data. That is, with the full set of 
data points, r]i at the respective temperature i = 1, 2, TV, where TV is the total number of supercooled liquid data points determined 
by locating T OI Error 2 = (1/iV) [log yjTj) — logr^] 2 , where r](T) is the fitting function at temperature T with parameters determined 

by minimizing the mean square deviation of the fitting function from experiment in the indicated range of experimental data. 



TABLE II: Comparison between double exponential and parabolic fit for supercooled liquid B2O3 





Double Exponential" 


Parabolic" 


Fitted temperature range 


log^^'/Poise 


K 


C 


Error 6 


J/K" 1 


T /K 


log 770/Poise 


Error 6 


533 K - 970 K 


0.524 


511 


1.41 • 10 3 


0.13 


3.51 • 10 3 


1.07 • 10 3 


2.96 


0.087 


675 K - 820 K 


1.61 


169 


2.01 • 10 3 


0.46 


3.46 ■ 10 3 


1.05- 10 3 


3.14 


0.20 


533 K - 675 K 


-0.940 


974 


1.12 • 10 3 


0.25 


3.74 • 10 3 


1.01 • 10 3 


3.14 


0.14 



Listed fitting parameters are determined by minimizing the mean square deviation between the fitting function and experiment for the 
specified range of data. 
6 As defined in Table [j] 



TABLE III: Comparison of Fractional Exponential and Parabolic fit for PPG 





Fractional Exponential" 


Parabolic" 


Fitted temperature range 


logr c (oo) / s 


X 


T c /K 


Error 6 


J/K- 1 


T /K 


lOgTo/s 


Error 6 


200 K - 240 K 


-5.49 


72.4 


244 


0.068 


2.26 ■ 10 3 


264 


-6.12 


0.042 


222 K - 240 K 


-6.66 


41.8 


265 


0.33 


2.27- 10 3 


265 


-6.25 


0.13 


200 K - 218 K 


-3.73 


116 


227 


0.11 c 


2.29 ■ 10 3 


265 


-6.39 


0.088 c ' d 



a Listed fitting parameters are determined by minimizing the mean square deviation between the fitting function and experiment for the 
specified range of data. 
6 As defined in Table [i] 

c When the fit range is restricted to the lowest temperatures, T c is lower than the lowest data point available for PPG. Therefore, the 
Error is calculated only for T < 227 K for both the fractional exponential and the parabolic fit for sake of comparison. This requires that 
the seven highest temperature points be excluded from consideration. If those points are included, the Error for the fractional exponent 
over the entire range considered would become undefined. 

d Over the entire range for which T <T = 264 K, the Error for the the parabolic fit using the parameters for this range becomes 0.13. 



Nagel 22] applied a generalization of Eq. [T] in this way 
to collapse seemingly disparate data on a soft-sphere sys- 
tem [53]. We further illustrate the organizational power 
in this paper by considering the role of attractive forces in 
supercooled transport. Such forces can produce large ef- 
fects at supercooled conditions [21] • We show that these 
effects reflect modest trends in the parameters J and T . 



II. COMPARATIVE STUDIES 
A. VFT form compared with parabolic form 

The VFT formula for transport property r is 

log(r/4^) = - % , (2) 

1 — J VFT 



where t V ft! A, and Tvft are fitting parameters. It is 
possibly the most common of all expressions used in glass 
physics, often referred to as a law [25]. This formula is 
often invoked to attribute slowing dynamics in a super- 
cooled liquids to a thermodynamic transition at the tem- 
perature Tvft- Adam-Gibbs theory and random- 
flrst-order theory |12] connect this transition to an en- 
tropy crisis, where Tvft is a recorder of a Kauzmann 
temperature. The latter, Xk, is a temperature at which 
the entropy of the supercooled liquid equals that of the 
ordered solid. (A formula analogous to Eq. [2] is an exact 
result for a class of stochastic models [26] .) 

It is often asserted that Tvft ~ Ik is a good approx- 
imation |12| . but it has been observed that experimen- 
tal evidence is mixed [27]. To analyze the issue quan- 
titatively, note that both Tvft and Tk are constrained 
to lie below the glass transition temperature, T g , and 
the relative difference from that reference temperature 
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FIG. 1: Examples of using Eq. [T] (a) or Eq. [5] (b) to fit 

viscosity, r), of a supercooled liquid as a function of recipro- 
cal temperature, 1/T. The circles are experimental data [29] . 
Three fits are shown for both equations. For one fit, parame- 
ters are determined by minimizing the mean square deviation 
between functional form and experiment for the full range of 
supercooled data, for the other two, parameters are found by 
minimizing the mean square deviation between the functional 
form and a subset of that data, the subset being either the 
higher temperature range of data or the lower temperature 
range of data. See Table [I] for specified ranges and param- 
eters. The arrow indicates the value of 1/T , marking the 
crossover between normal and supercooled liquid behaviors. 
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FIG. 2: Examples of using Eq. [T] (a) or Eq. [3] (b) to fit 
viscosity, r), of a supercooled liquid as a function of recipro- 
cal temperature, 1/T. The circles are experimental data [33] . 
Three fits are shown for both equations. For one fit, parame- 
ters are determined by minimizing the mean square deviation 
between functional form and experiment for the full range of 
supercooled data, for the other two, parameters are found by 
minimizing the mean square deviation between the functional 
form and a subset of that data, the subset being either the 
higher temperature range of data or the lower temperature 
range of data. See Table [TT] for specified ranges and param- 
eters. The arrow indicates the value of 1/T , marking the 
crossover between normal and supercooled liquid behaviors. 



is the pertinent measures. According to the tabulated 
temperatures for 33 different liquids |28| . the quanti- 
ties (Tg/TvFT — 1) and (T g /Tk — 1) are typically a few 
tenths. Their ratio, R = {T g /T VFT - 1) / (T g /T K - 1), 
often differs significantly from unity. In particular, 
averaging over the 33 different liquids yields a mean 
value of (R) = 1.15 and a root-mean-square deviation 



((R- (R)f) 1/2 = 0.613. In other words, it is just as 
likely that Tvft and Tj< will differ significantly as it is 
that they will be similar. 

Independent of an alleged connection to thermody- 
namics, the utility of the VFT formula can be tested 
by examining whether its parameters are material prop- 
erties. We do so for one typical fragile liquid, o-terphcnyl 
(OTP). What we illustrate with this system is represen- 
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tative of what we generally find for a wide selection of 
supercooled liquids. 

The onset temperature for OTP is about 350 K, and 
the glass transition temperature is 239 K. Data exists 
over this entire range [33]. Our earlier paper [5] con- 
sidered the lower half of this range, but here we con- 
sider the full supercooled range for the purpose of ex- 
ploring sensitivity of parameters. (We omit the normal 
liquid range T > T a , and only remark that the defi- 
ciencies of VFT we are about to demonstrate are more 
severe when the temperature range is extended to in- 
clude the normal liquid regime.) With this range of 
data, we have fit the parameters for the VFT formula 
and those for the parabolic formula in three different 
ways: by including data from only the range of higher 
supercooled temperatures (267 K < T < 350 K), by in- 
cluding data from only the range of lower supercooled 
temperatures (239 K < T < 267 K), and by including 
data from the entire range of supercooled temperatures 
(239 K < T < 350 K). The fitting parameters obtained 
in these ways are given in Table [I] 

The entries to this table show a significant difference 
between parameters found from fitting the full range of 
data and those found from fitting only the higher tem- 
perature data. The VFT parameters found from fitting 
only the lower temperature data agree well with those 
found from fitting all supercooled data at the expense of 
a larger fitting error. For these cases, the error in repro- 
ducing the full range of supercooled liquid data is rea- 
sonably small, but at the expense of an unphysical value 
for the reference viscosity, i]^^, • Figure 111 shows poor 
agreement between experiment and the VFT formula for 
high temperature data when only low temperature data 
is included in fitting (and vice versa). In the latter case, 
when only the higher temperature range of supercooled 
data is used to fit VFT parameters, the value of Tvft is 
close to T g , causing the huge error reported in the second 
row of Table [I] In contrast, the parameters and quality of 
fits found with the parabolic form change little between 
the full set of data or either of the subsets of data. 

Thus, parameters of the parabolic form appear to rep- 
resent properties of the system, and those determined 
over one range of temperatures can be used to reliably 
predict properties over another range of temperatures. 
The parameters of the VFT form, however, depend upon 
both the properties of the system and the range of data 
considered, and as such the VFT form cannot predict 
properties outside the range over which it has been fit. 



B. Double-exponential form compared with 
parabolic form 

The double-exponential formula for transport property 
r is 
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FIG. 3: Examples of using Eq. I] (a) or Eq. |4](b) to fit relax- 
ation time, r, of a supercooled liquid as a function of recipro- 
cal temperature, 1/T. The circles are experimental data [34] . 
Three fits are shown for both equations. For one fit, parame- 
ters are determined by minimizing the mean square deviation 
between functional form and experiment for the full range of 
supercooled data, for the other two, parameters are found by 
minimizing the mean square deviation between the functional 
form and a subset of that data, the subset being either the 
higher temperature range of data or the lower temperature 
range of data. See Table |TTT| for specified ranges and param- 
eters. The arrow indicates the value of 1/T , marking the 
crossover between normal and supercooled liquid behaviors. 
The leftmost endpoints of the fit lines in figure (b) indicate 
the upper-temperature end point for applying Eq. [4] This 
temperature changes depending upon the range of data con- 
sidered, and data for T > T c must be excluded from fits using 



where 



(3) 



dx , K, and C are fitting parameters. A special 
case, where K and C are of the same order, is the behav- 
ior of two-spin (or two-particle) facilitated lattice models 
[30] . Over the lower half of the temperature range be- 
tween T and T g , the double exponential has been used to 



5 



successfully collapse transport data [3T]. Mauro and co- 
workers [19] and others [32] have applied it to a broader 
range and report that Eq. [3] is superior to both the VFT 
form, Eq. [2] and the parabolic form, Eq. [T] In fact, 
the double-exponential form suffers from the same mal- 
ady as the VFT form when applied to fit data over the 
full range of supercooled temperatures. As a data set is 
enlarged, its parameters fail to converge, implying these 
parameters are not material properties and the formula 
cannot be used to predict data not yet measured. The 
problem is not as severe as it is for the VFT form, owing 
to the fact that the double exponential is not singular 
while the VFT expression is singular. But the deficiency 
of the double-exponential form is nonetheless significant, 
as illustrated in Figure [2] and Table [IT] 

Our illustration considers the inorganic glass forming 
liquid B2O3 [33]. The behaviors found for this system are 
typical of what we find for several other systems. For this 
particular system, the onset temperature for this liquid 
is close to 1000K, and the glass transition temperature 
for this liquid is about 500K. Fitting data over this entire 
range, the double-exponential form proves reasonably ac- 
curate, but its parameters change markedly as the range 
of fitted data changes. Fitting only higher temperature 
data, 675 K sC T ^ 820 K , yields a function that inac- 
curately describes the lower temperature data, and fit- 
ting only lower temperature data, 533 K ^ T ^ 675 K, 
yields a function that inaccurately describes the higher 
temperature data. In contrast, the parameters and ex- 
cellent quality fits of the parabolic form change little as 
the range of fitted supercooled data changes. 

Thus, for the purpose of employing a particular func- 
tional form to predict low temperature properties from 
measured properties at higher supercooled temperatures 
(or vice versa), Eq. [3] is superior to the VFT expression, 
but this double exponential form advocated by Mauro 
and coworkers |19j is inferior to the parabolic form, Eq. 

ED 

C. Fractional-exponent form compared with 
parabolic form 

In both the VFT and double-exponential forms, tem- 
perature variations are more rapid than in the parabolic 
form. What if temperature variation is taken to be 
less rapid? An alternative of this type is the specific 
fractional-exponent form suggested by numerical solu- 
tions to Saltzmann and Schweizer's theory for structural 
relaxation in polymer melts [20) . 

log(r/r c )=X(r c /T-l) 1 - 57 . (4) 

We illustrate the performance of this expression in Figure 
[3] and Table |III| We do so with data for the glass form- 
ing polymer melt polypropylene glycol (PPG) [34]. The 
onset temperature for this liquid is about 264 K, and its 
glass temperature is about 199 K [5]. The data covers 
most of this supercooled region, but not all the way up 



to the onset. As in the previous two case studies, we con- 
sider three ranges of the existing data: the entire data 
set, which extends up to about 10% of the onset tem- 
perature, T Q , a lower temperature subset of that data, 
and a higher temperature subset of that data. As with 
the VFT and double-exponential forms, we find that the 
fitting parameters for this fractional-exponential form de- 
pend upon the range of data considered. Our illustration 
of this fact is typical of what we find when treating other 
materials with the same analysis. 

Due to the nature of the fractional exponent, this func- 
tional form can be used to fit data for T > T c . If the ex- 
ponent were 2, rather than 1.57, this temperature would 
be the onset temperature. When using an exponent of 
1.57 with subsets of the data where temperatures are all 
much smaller than T Q often one obtains T c < T , so that 
less than the complete set of supercooled data can be 
covered. Further, unlike the parabolic form with param- 
eters that vary little with changing data sets, the best fits 
of Eq. [4] produce parameters T c and X that vary widely 
from one subset of data to another. Moreover, the func- 
tion obtained fitting parameters to the higher tempera- 
ture range of data, 222 K ^ T ^ 240 K, provides an in- 
accurate representation of the data at the lower tempera- 
tures, and the function obtained fitting parameters to the 
lower temperature range of data, 200 K ^ T ^ 218 K, 
provides an inaccurate representation of the data at the 
higher temperatures. In contrast, the excellent quality 
of fits obtained with the parabolic form, Eq. 1, change 
little as the range of fitted data change. 

Thus, the fractional exponent 1.57 proves to be less 
satisfactory than that of the parabolic form, Eq. [T] 



III. ROLE OF ATTRACTIVE FORCES 

The prior section presents evidence in support of Eq[T] 
as a universal form for transport properties of fragile 
supercooled liquids. In this section, we show that by 
accepting the validity of this form we sensibly organize 
recent simulation results that might appear puzzling in 
the absence of this organization. In particular, Berthier 
and Tarjus[2"3] have shown that at some supercooled tem- 
peratures and densities relaxation in the Kob-Andersen 
Lennard- Jones (LJ) mixture[35 is orders of magnitude 
slower than that in the corresponding Weeks-Chandler- 
Andersen[36] (WCA) mixture. The difference between 
the WCA and the LJ potential reflects the significance of 
attractive forces, which are present in the LJ mixture and 
absent in the WCA mixture. Does the Berthier- Tarj us 
finding reveal a new mechanism for glassy physics, one 
that does not follow from constraints and local rigidity 
imposed by repulsive forces? 

The effects of attractive interactions uncovered by 
Berthier and Tarjus appear in dynamics, but not struc- 
ture. They show that the pair distribution functions of 
the two mixtures differ very little. Dynamics is associ- 
ated with fluctuations away from mean local structure. 



6 




FIG. 4: (a) parabolic collapse for WCA & LJ binary mixture simulation data for various densities, p from [21]. Parameters 
of fit are described in Eq. [T] Fitting parameter trends using Eq. [I] for the Kob- Andersen LJ and corresponding WCA binary 
mixtures from [24] for various net particle densities, p. (b) The inverse onset temperature, 1/T , as a function of density p. 
(c) Logarithm of relaxation time at the onset temperature, logr , as a function of p. (d) Transport energy parameter, J, as a 
function of p. Dotted line in (a) is the universal parabolic form. Lines connecting points in (b), (c) and (d) are guides to the 
eye. Error estimates are the size of the symbols. The unit of time is At = (m<7 AA /48eAA) 1//2 . 



Repulsive forces dominate the most likely equilibrium ar- 
rangements. They also constrain motions of particles, so 
that any spatial reorganization is rare and requires the 
coordinated displacements of several particles. As such, 
potential energy barriers to reorganization will have con- 
tributions from the interactions between several pairs of 
particles. The sums of small contributions from many 
particles can become significant. In other words, it can 
be due to the local rigidity from packing forces that small 
effects from attractive forces can become notable. 

To pursue this idea, we have examined the effects of at- 
tractive forces on supercooled fluid transport through the 
behaviors of J and T Q because these parameters charac- 
terize the energetics of these collective displacements. [T3] 
We have fit the simulation data[23 reported to determine 
the values of these parameters as functions of the liquid- 
mixture density. These results are shown in [4] along with 
the reference time scale r Q . While the values of J, T , 
and log t do vary slightly between the two types of mix- 



tures at the same density, their overall trends as functions 
of density are the same. Both liquid mixtures therefore 
seem to behave as similar supercooled liquids. 

Berthier and Tarjus remark on the "enormous" factor 
of 10 3 ' 2 by which r for the LJ mixture differs from that 
for the WCA mixture at pa AA = 1.2 and ^bT/eaa = 
0.41. This effect reflects modest differences in J/caa 
and /cbTo/caa, which are 1.8 and 0.73 for the LJ mix- 
ture, and 0.92 and 0.60 for the WCA mixture, respec- 
tively. The onsets to supercooled behavior thus appear 
at similar temperatures, and the activation energies for 
dynamics differ by less than one intermolecular attractive 
energy. Moreover, as the density increases the character- 
istic parameters for the two systems seem to converge to 
the same respective values. These findings are consistent 
with empirically established fact that the basic princi- 
ples of supercooled liquids are well captured by models 
without attractive intermolecular forces [37] [35] . 

Thus, while enhanced by the results of Berthier and 
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1000 K/T 

FIG. 5: An example of fitting parameter convergence for 
B2O3. Red points indicate the minimal data set required to 
converge fitting parameters when beginning from the low tem- 
perature data and extrapolating downward. Arrows indicate 
the location of the onset of fragile behavior. 

Tarjus, the picture of glassy physics as a class of phe- 
nomena caused by local constraints and rigidity remains 
unchanged. 
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Appendix: How to fit data to the parabolic form. 

Here, we summarize the procedure we employ to fit 
data with the parabolic form. 

1. Relaxation time or viscosity data spanning several 
orders of magnitude are required. 

2. Examine the data to attempt to identify normal liquid 
and fragile regimes, and thereby obtain an approximate 
location of the onset lying between the two. 

3. Fit data spanning several orders of magnitude start- 
ing with the lowest temperatures available, using a least 
squares analysis to obtain a first estimate of J, T and 
logT . 

4. Use this estimate for T a as an indication of the highest 



temperature point for which to fit the data and continue 
to add higher temperature points to your fit until the 
values of J, T and logr Q have converged. (An illustra- 
tion of the minimal amount of data needed to fit viscosity 
data for B2O3 [33J is shown in Figure [5j An illustration 
of the parameter convergence for the B2O3 fit is shown 
in Figure [6j) 

As data points are added after convergence, small 
systematic changes in parameter values can occur. 
Generally, J and T Q tend to rise together where as log r] 
decreases and vice versa. These variations reflect the 
presence of local minima in fits of nearly equal standard 
deviation that swap from being local to global, and 
their sizes appear to be negligible in comparison with 
statistical errors in experimental data. 

Caveats: 

A. While in our B2O3 example we have chosen to start 
at the lowest temperature point and add higher and 
higher temperature points until the parameters converge, 
it should be noted that often the lowest temperature 
points are the hardest to measure and therefore some- 
times data in this region is not always as reliable as data 
taken from a more moderate, yet still supercooled regime. 
It may be useful to check if a fit through such a moderate 
range properly predicts the lowest temperature data as 
an estimate for the quality of the lowest temperate data. 

B. The idea that fragile glass former transport data may, 
at low temperatures, crossover again from fragile to Ar- 
rhenius has been proposed |14j and recently considered 
further [39]. At this moment, there is limited evidence to 
support or refute this claim. Certain systems seems to 
exhibit this phenomenon [39-41 . However, the systems 
considered are either confined to nano-scale pores, which 
may or may not reflect the behavior of macroscopic ma- 
terials, or more recent experiments contradict those that 
exhibit a crossover 42J. In the future, it may be neces- 
sary to reconsider this idea if lower temperature data can 
be collected for experimental systems. In this case, for 
the parabolic fitting, extra caution should be taken with 
the data and a low temperature fragile-strong crossover 
point T x should be identified. Under these circumstances, 
fit the lowest temperature points with: 

log(r/. x ) = 

From this, choose a reasonable low-temperature end to 
the fragile fitting regime, T x , where the Arrhenius behav- 
ior at low temperatures diverges from the experimental 
or numerical data. 
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FIG. 6: An example of fitting parameter convergence for B2O3. In the fragile regime, 47 data points ranging continuously over 
11 orders of magnitude in r\ are available for fitting. Including the liquid regime, there are 62 data points in this set spanning 
13 orders of magnitude in 77. Initially, only the lowest temperature points are considered and subsequent fits with more points 
come from adding higher temperature points to the considered set. (a) shows the convergence of T a as a function of number 
of points used for fitting, (b) shows the convergence of J as a function of number of points used for fitting, (c) shows the 
convergence of r\ (analogous to r ) as a function of number of points used for fitting. For this system, T — 1066 K, J/T = 
3.3, and log?7 /Poise = 2.9. 



[1] M. Ediger, C. Angell, and S. Nagel, J. Phys. Chem. 100, 
13200 (1996). 

[2] K. Binder and W. Kob, Glassy Materials and Disordered 
Solids (World Scientific, New Jersey, 2005). 

[3] C. Angell, K. Ngai, G. McKenna, P. McMillan, and S. 
Martin, J. Appl. Phys. 88, 3113-3157 (2000). 

[4] C. Angell, Science 267, 1924-1935 (1995). 

[5] Y. Elmatad, J. P. Garrahan, D. Chandler, J. Phys. 
Chem. B 113, 5563 (2009). 

[6] J. Jonas, Science 216, 1179 (1982). 

[7] D. Chandler, J. D. Weeks, H. C. Andersen, Science 220, 
787 (1983). 

[8] H. Tyrrell, K. Harris. Diffusion in Liquids: A theoretical 
and experimental study (Butterworth Publishers, Stone- 
ham, 1984). 

[9] D. Chandler and J. P. Garrahan, Annu. Rev. Phys. Chem. 

61, 191-217 (2010). 
[10] G. Biroli and J. -P. Bouchaud, The Random First-Order 

Transition Theory of Glasses: a critical assessment, 
|arXiv:0912.2542| (2010). 
[11] G. Adam and J. Gibbs, J. Chem. Phys. 43, 139 (1965). 
[12] V. Lubchenko and P. Wolynes. Annu. Rev. Phys. Chem 

58, 235-266 (2007). 
[13] P. Sollich and M.R. Evans, Phys. Rev. Lett. 83, 3238 

(1999). 

[14] J. P. Garrahan, D. Chandler, Proc. Nat. Acad. Sci. 100, 
9710 (2003). 

[15] F. Ritort and P. Sollich, Adv. Phys., 52, 219 (2003). 
[16] L. Berthier and J. P. Garrahan, J. Chem. Phys. B 109, 
3578 (2005). 

[17] D. Ashton, L.O. Hedges and J. P. Garrahan, J. Stat. 

Mech. P12010 (2005). 
[18] G. Tammann, J. Soc. Glass Technol. 9, 166 (1925). 
[19] J. Mauro, et al, Proc. Nat. Acad. Sci. 106, 19780 (2009). 
[20] E. Saltzman and K. Schweizer, J. Chem. Phys. 121, 2001 



(2004). 

[21] D. Kivelson, S. A. Kivelson, X. L. Zhao, Z. Nussimov abd 

G. Tarjus, Physica A 219, 27 38 (1995). 
[22] N. Xu, A. Liu and S. Nagel, Phys. Rev. Lett. 103, 245701 

(2009). 

[23] L. Berthier, T. A. Witten. Phys. Rev. E 80, 021502 
(2009). L. Berthier, T. A. Witten Europhys. Lett. 86, 
10001 (2009). 

[24] L. Berthier, G. Tarjus. Phys. Rev. Lett. 103, 170601 
(2009). 

[25] A recent example: M. Wyart Phys. Rev. Lett. 104, 
095901 (2010). 

[26] C. Toninelli, G. Biroli, and D.S. Fisher. Phys. Rev. Lett. 

96 035702 (2006). 
[27] H. Tanaka. Phys. Rev. Lett. 90, 055701 (2003). 
[28] C. A. Angell, J. Res. Natl Inst. Stand. Technol. 102, 171 

(1997). 

[29] D. Plazek, C. Bero, I.-C. Chay, J. Non-Cryst. Solids. 

172-174, 181 (1994). 
[30] C. Toninelli, G. Biroli, and D.S. Fisher. J. Stat. Phys. 

120 167-238 (2005). 
[31] E. Rossler, K.-U. Hess, V. Novikov, J. Non-Cryst. Solids. 

223, 207 (1998). 
[32] P. Lunkenheimer, S. Kastner, M. Kohler, and A. Loidl, 

Phys. Rev. £81, 051504 (2010). 
[33] H. Tweer, J. H. Simmons, P. B. Macedo, J. Chem. Phys. 

54, 1952 (1971). 
[34] N. B. Olsen, T. Christensen, and J. C. Dyre, Phys. Rev. 

Lett. 86, 1271-1274 (2001). 
[35] W. Kob and H. Andersen, Phys. Rev. Lett. 73, 1376 

(1994). 

[36] J. Weeks, D. Chandler, and H. Andersen, J. Chem. Phys. 
54, 5237 (1971). In the context of glassy physics models, 
WCA mixtures have been first studied by: L. Maibaum. 
Fluctuations and dynamics in the liquid state. Doctoral 



9 



Thesis, UC Berkeley, Berkeley, 2005; L. O. Hedges, L. 

Maibaum, D. Chandler, J. P. Garrahan, J. Chem. Phys. 

127, 211101 (2007). 
[37] N. Gnan, T. B. Schr0der, U. R. Pedersen, N. P. Bailey, 

J. C Dyre, J. Chem. Phys. 131, 234503 (2009). 
[38] U. R. Pederse n, T. B. Schr0der, J. C. Dyre, 

larXiv:1004.53T6l 
[39] S.-H. Chen, Y. Zhang, M. Lagi, S. Chong, P. Baglioni 

andF. Mallamace, J. Phys: Condens. Matter, 21, 504102 

(2009). 



[40] S.-H. Chen, et al, Proc. Nat. Acad. Set. 103, 9012 (2006). 
[41] W. T. Laughlin, D. R. Uhlmann, J. Phys. Chem. 76, 
2317 (1972). 

[42] See for example: P. Dixon, N. Menon, S. Nagel, Phys. 
Rev. E 50, 1717 (1994); C. Gainaru, et al, Phys. Rev. 
B. 72, 174203 (2005); R. Richert, C. A. Angell, J. Chem. 
Phys. 108, 9016 (1998); R. Richert, J. Chem. Phys. 123, 
154502 (2005). 



